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Increasingly emerging technologies in agriculture such as computer vision, 
artificial intelligence technology, not only make it possible to increase 
production. To minimize the negative impact on climate and the 
environment but also to conserve resources. A key task of these technologies 
is to monitor the growth of plants online with a high accuracy rate and in 
non-destructive manners. It is known that leaf area (LA) is one of the most 
important growth indexes in plant growth monitoring system. Unfortunately, 
to estimate the LA in natural outdoor scenes (the presence of occlusion or 
overlap area) with a high accuracy rate is not easy and it still remains a big 
challenge in eco-physiological studies. In this paper, two accurate and non- 
destructive approaches for estimating the LA were proposed with top-view 
and side-view images, respectively. The proposed approaches successfully 
extract the skeleton of cucumber plants in red, green, and blue (RGB) 


images and estimate the LA of cucumber plants with high precision. The 
results were validated by comparing with manual measurements. The 
experimental results of our proposed algorithms achieve 97.64% accuracy in 
leaf segmentation, and the relative error in LA estimation varies from 3.76% 
to 13.00%, which could meet the requirements of plant growth monitoring 
systems. 
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1. INTRODUCTION 

Modern agriculture with high technology has been played a crucial role in feeding the growing 
population [1]. By efficiently utilizing the limited resources, using modern technologies in farming, farmers 
in modern agriculture enhance both quantity and quality of agricultural products [2]-[4]. Plant growth 
monitoring based on machine vision is one of such modern technologies [5]. This monitoring approach using 
digital cameras and suitable sensors for image acquisition. The growth analysis of plant phenotypes in a non- 
destructive manner is conducted automatically [6]-[8]. However, some plant features, such as stem diameter, 
leaf number, and leaf area are not measured easily by such non-contact measurement methods [9]. Therefore, 
improving the quality of plant features estimation is one of the primary concerns of present-day plant 
management [10]-[15]. 

Leaf is one of the most important photosynthetic organs of a plant [16]. Hence, leaf area is an 
important factor in measuring the growth and productivity of plants. Consequently, the precise determination 
of leaf area is basic in many plant physiological studies [17]-[20]. In fact, the leaf area is traditionally 
measured by regression equation methods [21], grid counting methods [22], gravimetric method [23], and 
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some leaf area meters (planimeter) [24], [25]. However, most of these methods are destructive (they require 
the excision of leaves, which may damage the growth of plants), and provide off-line performance analysis. 
Though these traditional methods have high accuracy, they require expensive equipments and time- 
consuming experiments. Fortunately, some nondestructive methods for LA estimation [26]-[29] have been 
being developed rapidly to minimize the disadvantages of these destructive methods. It is proven that the 
nondestructive methods achieve better performances, such as faster, cheaper, and provide real-time 
experimental analysis than the destructive methods. 

Nowadays, nondestructive measurement methods based on digital image analysis are widely used in 
many growth models. The authors in [27] have focused on measuring cucumber parameters based on 
machine vision technology. It is a real non-destructive online measurement method that used digital cameras 
to capture the canopies of cucumber plants. The image parameters of cucumber plants are combined with 
manual measurements to construct some inversion models for cucumber growth parameters. The results show 
that with high coefficients of determination between manual measurements and inversion values of cucumber 
parameters, these inversion models are able to be a non-destructive online method for measuring cucumber 
parameters. 

Wang et al. [26] developed a new image segmentation method based on the Chan-Vese (C-V) model 
and the Sobel operator. In order to improve the segmentation precision of overlapping cucumber leaves, the 
procedure of the image segmentation method is described as follows: First, the background is mitigated based 
on a threshold with respect to the relative levels of green in the image. Using the C-V model, the contour of 
the target leaf in the upper layer is extracted. Then, the leaf edges of the target leaf are detected by a gradient 
operator with eight directions. Finally, the combination of the results obtained by the C-V model and the 
Sobel operator is used to extract the overlapping target leaf. The experimental results show that with a low 
level of the mean error rate (0.0428), the proposed method is successful in extracting the target leaf from 
complicated background images of cucumber plants. 

Ibrahim et al. [30], two image segmentation methods such as k-mean clustering and region growing 
have been used for automatic flood detection. The simulation results reveal that the region growing using 
gray images achieved higher segmentation accuracy than that of k-mean clustering method. In [31], a system 
for detection and classification of various pest attacks and infection on plants is introduced. For segmentation 
of preprocessed sample, an image segmentation method based on hybrid watershed segmentation with extended 
K-Means clustering is proposed. The experimental results show that the proposed approach achieved higher 
accuracy rates than some existing approaches in disease/pest detection. More recently, Liu et al. [32] proposed 
an instance segmentation method, which used an improved version of mask region-based convolutional 
neural network (Mask RCNN) to detect cucumber fruits in pixel level. The experimental results reveal that 
the proposed improved Mask RCNN achieves high accuracy in detecting cucumber fruits. The precision and 
recall rates can reach 90.68% and 88.29%, respectively. All studies mentioned above show that LA of 
cucumber plants can be determined by such nondestructive methods with high accuracy. This important LA 
information helps growers having proper time irrigation, cultivation, therefore, get better yields. In this paper, 
we have focused on improving the accuracy of cucumber leaf area estimation, which could provide a 
valuable resource for managing the cucumber plant. 


2. RESEARCH METHOD 

According to Figure 1, the flowing chart of cucumber leaf area estimation is used to improve the 
research and experiment comparison. The flowchart illustrates the various steps of a process that a captured 
image is used for LA estimation. There are two different estimation techniques for top-view and side-view 
images. Each step of the process is described in more detail in following subsections. 


2.1. Background elimination 

To improve the segmentation accuracy of overlapping target leaves, some approaches [26], [33] 
have focused on removing background (the non-leaf areas such as soil, stems, and petals). On considering the 
difference of colour between background areas and foreground areas (leaf areas) [33], the authors in [26] 
used the levels of the green channel in the RGB colour space to remove all pixels in non-leaf areas. In that 
study, the relationship of pixel colour channels for a green leaf is determined in (1). 


(Gratue = Byatue > 01) N (Gratue = Ryatue > 02) (1) 


where Rygive, Byaiue, aNd Gyaine are value of red, blue, and green component, respectively. 04, Ozare the 
Ryatue and Byatue control parameters, respectively. These parameters vary in different species of plant, e.g., 
they are chosen 6, = 10, and 0, = 15 for leaf segmentation of cucumber leaves. In this way, we could 
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remove all non-leaf areas by just keeping pixels whose relationship of colour channel is green. Unfortunately, 
if the difference between the background colour and the foreground colour is not apparent or the colours of 
leaves change in the autumn, this approach might not work properly. 


Image preprocessing 


Background elimination 


Top-view images 


Estimated process 


a 


Validation process 


Figure 1. The flowing chart of cucumber leaf area estimation 


Fortunately, the Mahalanobis distance [34] is an efficient alternative approach for background image 
removal that mitigates these above-mentioned limitations. The basic idea of this image segmentation method 
is to compute the distance between colour intensities of pixels and their neighbouring pixels within an 
ellipsoidal region. One pixel becomes a member of a group if its Mahalanobis distance is smaller than a 
threshold. In this way, the Mahalanobis distance based method for image background removal is expected to 
improve the quality and reduce the false alarms [35]. 

In fact, many Mahalanobis distance based image background removal algorithms are proposed [36], 
[37]. In [35], for weed leaf image segmentation, each pixel of the image has been classified as a weed area 
pixel or not based on the Mahalanobis distance from the plant colour class and a colour threshold value. The 
experimental results show that 84% of the weed leaves were correctly segmented. It is clear that this 
segmentation accuracy is quite high and acceptable for weed segmentation with multiple leaves in an image. 
However, not many of these studies consider how to determine an optimal threshold of each sub-region for 
leaf segmentation. It might increase the accuracy of image segmentation [38]. To do this, in this paper, we 
propose a Modified Mahalanobis Distance based method (MMD) for background image removal. 

Without loss of generality, let denote the color value of a pixel in a color image. We assume that the 
color of each pixel is the normal distribution. 

In order to build the foreground model, some initial frames that have only foreground (the leaf 
image of considered plants) are analyzed. Let denote the number of such initial frames. The mean of the 
colour at every pixel is determined in the (2) [39]. 


1 
HE Èk Xk (2) 
And the covariance is computed as follows: 


C= = Dee — Hk) (xk — My)" (3) 
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after making the background model, the Mahalanobis distances of every pixel are calculated using (4): 


d = | (x — u)T Ct (x — u) (4) 


Check if d > 09, where Oois predefined threshold, that pixel is classified as foreground, otherwise, it 
belongs to background. The problem here is how to determine the optimal threshold (65) for leaf 
segmentation of different plant?. Xia et al. [37] 65 = 10,40,7 respectively for pepper, old paprika, and young 
paprika after testing leaf images of these plants. Unfortunately, these leaf segmentation approaches based on 
a global threshold of the entire image may lead to low accuracy and a high rate of false-positive alarms. 
Therefore, in order to improve the accuracy and effectiveness in leaf segmentation, the threshold should 
change adaptively depending on properties of each sub-image. To do this, the authors in [40] implemented a 
new segmentation algorithm with the following steps. First, the standard deviation of the grey level of the 
whole image, which is named global standard deviation, is computed. Then the image is divided into four 
equal regions. Each sub-region is continuously split into four sub-regions if its standard deviation of the grey 
level (namely local standard deviation) is higher than the global standard deviation. This procedure is 
recursively continued until the local standard deviation is less than or equal to the global standard deviation 
of image grey level. The threshold of each sub-region is determined based on the local properties of the 
image by using IsoData algorithm [41]. 

In this study, we combine the Mahalanobis distance algorithm and optimizing threshold algorithm in 
[40] to propose a modified Mahalanobis distance-based method for background image removal. The details 
of our proposed MMD algorithm are described in algorithm 1. 


Algorithm 1: MMD algorithm 


Input: 

- Raw images. 

- A grey-level histogram hist. 

Output: 

- Optimal threshold of each sub-region 45 

- Segmentation images 

Initiation: 

Let X = {x1 X2... Xyxył} denote the data set of NxN pixels of source image file; Let L=256 
denote the number of grey levels in NxN image; n,is the number of occurrences of the grey 
levels; The local standard deviation g,=©. 

Procedure 

- MMD algorithm(L,n,, 69). 

- Calculate the global standard deviation oa: 


o= [ELSU — mye (5) 


Divide the original image D into four equal regions: 
D,,t = 1,2,...4;D = X Di; Di AD Vit) 


for (i = 1 to 4) do 
while 0; >ø do 
- Divide the current sub-region D; into four equal sub-regions Djj,j = 1,2,...,4 

- Calculate the local standard deviation of the current sub-region gj © Gj 

end while 

= Calculate the optimal threshold of current sub-region 09; = 0; 

- Calculate the Mahalanobis distance of investigated pixel (d;)by using (4) and 
compare its Mahalanobis distance with optimal threshold of sub-region 6 ;where 
the investigated pixel falls inside. 

if dj >0$; then the investigated pixel is classified as foreground 
else the investigated pixel belongs to background 
end if 
end for 
end procedure 


For evaluating the quality of our proposed image segmentation algorithm, manual segmentation 
images were used to compare with the results of this segmentation algorithm. The following performance 
indices have been used: 

a) There are four possible outcomes when classifying a pixel: true positive, false positive, true negative 
and false negative [42]. They are described more in detail in Figure 2. 
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- True positive (TP): Pixels of leaf areas are classified correctly as leaf pixels. 

False positive (FP): Pixels in non-leaf areas are misclassified as leaf pixels. 

- True negative (TN): Pixels in non-leaf areas are classified correctly as non-leaf pixels. 
False negative (FN): Pixels of leaf areas are misclassified as non-leaf pixels. 

b) The false negative rate (FNR) is computed as given in (6). 


FNR = (6) 
TP+FN 
c) The false positive rate (FPR) is computed as given in (7). 
FPR = (7) 
FP+TN 
d) The error rate (ER) is mathematically defined in (8). 
ER = FN+FP (8) 
TP+TN+FP+FN 


Ground truth Segmented 
image image 


22242424 


NY 
NI 
RNN 


Figure 2. Schematic illustrating the definition of four possible outcomes when classifying a pixel 


2.2. Measurement of leaf area 

Nondestructive leaf area estimation with a high accuracy rate and in a small running time is critical 
in many specific fields of plant growth modelling (e.g., plant physiological and ecological experiments [43]). 
Leaf area may be used for leaf area index calculation [9], [44] which is likely to model vegetation 
productivity and photosynthetic capacity [45], [46]. Moreover, LA is a very useful parameter in plant growth 
and morphological studies of the plants [47]. 

In this study, after removing the background, the leaf areas are measured by counting leaf pixels and 
comparing them with red calibration pixels. In this way, our proposed approach is able to estimate the leaf 
areas automatically while mitigating the need for manual operations for camera distance calculations. 
However, there are some overlap areas of the leaves, which may be considered as one layer of leaves. 
Additionally, the difference between green leaves and green objects may not be distinguished correctly if we 
only use top view images of the plant. Therefore, in our proposed algorithm, we use both top view and side 
view images of the plants to estimate the leaf area. 

The following steps are used to estimate the leaf area of the cucumber plant in the top-view image. 
Step 1: Convert the original image into a binary image. 

Step 2: Count the number of pixels in the calibration area (P,,) (The size of the calibration area given in this 
paper is 25cm?). 
Step 3: Count the number of pixels in the leaf image (P,) of the considered plant. 
Step 4: The leaf area can be obtained by (9). 
LA, =! (9) 
Pel 

The side-view images of plants are used to estimate the main stem length of plants and then to 
estimate the number of leaves. According to the number of leaves, we estimate the leaf area of the plant. The 
following steps are used to estimate the leaf area of cucumber plant in the side-view image: 

Step 1: Convert the original image into a binary image. 
Step 2: Skeleton extraction of cucumber plant by using Zhang-Suen algorithm [48]. The main stem of 
cucumber is described in white pixels. 


Digital image processing methods for estimating leaf area of cucumber plants ... (Uoc Quang Ngo) 


322 o ISSN: 2502-4752 


Step 3: Calculate the length of the main stem of cucumber plant: Scan the image from the bottom-left corner 
to the top-right corner of the image frame to catch and locate all white pixels. The length of the main 
stem of the cucumber plant is cumulative sum of Euclidean distances between these white pixels. 

Step 4: Count the number of nodes along the stem of the cucumber plant. And then estimate the number of 
leaves on the main stem. 

Step 5: According to the results given in [49] and our manual measurements given in Table 1, the estimated 
leaf area (LA,) is the cumulative sum of leaf area from the oldest leaf (the first leaf) to the youngest 
leaf on the cucumber stem. 


Table 1. Measured data 


Average of leaf area (cm’) 


par ot L 
number 1st 2nd 3rd 4th Sth 6th 7th 8th 9th 10th llth 12th LAm Cm 
10 2 07.25 100.07 207.32 12.35 
12 25 08.44 100.70 46.25 255.39 18.21 
14 3.62 10.92 101.20 92.48 66.51 371.10 23.89 
16 4.45 11.97 101.93 94.17 90.05 64.82 462.94 30.21 
18 5.61 13.72 103.21 95.22 89.91 74.38 58.39 534.83 35.72 
20 6.45 15.74 104.91 96.68 92.35 75.36 80.38 48.14 613.56 41.55 
22 7.65 21.26 106.30 98.10 92.96 77.38 97.13 90.99 78.81 762.92 47.42 
24 8.8 27.95 106.73 98.41 93.65 78.03 110.26 93.12 108.91 90.42 907.50 53.71 
26 9.75 33.82 122.70 112.94 110.99 88.85 117.83 103.46 118.24 121.47 83.98 1114.27 59.29 
28 10.92 37.68 130.18 136.44 118.38 103.12 119.98 108.30 126.41 128.70 116.52 102.66 1328.38 64.35 


30 11.95 41.61 135.96 147.32 131.55 129.33 123.94 123.98 135.20 135.86 127.79 130.02 115.4 1578.01 70.65 


2.3. Experiments 

Image acquisition system: The cucumber plants are labelled daily with QR codes which contain 
species, day after transplanting of cucumber. As depicted in Figure 3, our machine vision system consists of 
two cameras (Nikon camera D3500 DX), a personal computer (Core 17-7700 CPU, OS: MS Windows 10) 
and our image processing program is developed in Python software. The top-view image acquisition sytem in 
Figure 3(a) is a checkerboard target with known dimension is used as calibration area. The other calbiration 
area is also used the side-view image acquisition sytem in Figure 3(b). As depicted in Figures 4(a) and (b), all 
red squares are used as referenced area for leaf area estimation. 

Experimental design; The experiments are conducted from March 1 to April 30, 2021 in an East- 
West, polyethylene-covered greenhouse at the Vietnam National University of Agriculture 
(lat.21.0079764’N, long. 105.9335549E), Vietnam. The experimental greenhouse had two plots and 200 
cucumber plants (Cucumis sativus L. cv. Deltastar). The distance between two adjacent plants is 30 cm. Four 
plants in one plot are selected randomly for vegetation measurements. Images were taken on sunny days in 
direct sunlight. The image capture time was once at 11h:30 AM in every two days. The experimental images 
are shown in Figure 4 with top-view image in Figure 4(a) and side-view image in Figure 4(b). 


l 2000 mm 


1000 mm 


(b) 


Figure 3. Block diagram of image acquisition system: (a) top-view, and (b) side-view 


Indonesian J Elec Eng & Comp Sci, Vol. 25, No. 1, January 2022: 317-328 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 o 323 


Figure 4. Samples of cucumber plant images: (a) top-view, and (b) side-view 


In this paper, we used the manual measurements which serve as references for evaluating the 
performance of our proposed algorithms. The manual measurements included leaf number, the length of the 
main stem (Lm), the width and length of each leaf. According to the results given in [49], in this paper the 
measured value of leaf area (LAm )was carried out as given in (10). 


L; = —29.96 + 0.849 * lws * lis 
LAm = Xs L; (10) 


Where lws and l are measured width and length of leaf at position ston cucumber plant, 
respectively. 80 cucumber plants were randomly selected for manual measurement. Leaves at different 
positions (s from 1st (the oldest leaf) to 12th leaf) are measured 50 times and the average values of these 
parameters are given in Table 1. 


3. RESULTS AND DISCUSSION 
3.1. Segmentation results 

Eighty cucumber plants were used for leaf segmentation. In Figure 5, the segmentation results of 
top-view images are shown. It can be seen that the original image in Figure 5(a), is divided into some sub- 
regions. In each sub-region, an optimal threshold is chosen for comparison with the mahalobis distance of all 
pixels falling inside this sub-region. 

In Figure 5, there are 3 separate layers with different colour ranges. The first layer is the leaf area of 
cucumber plants in Figure 5(b). In this layer, all pixels were considered as part of cucumber leaves. 
The second layer denotes the calibration area in Figure 5(c), which is used for leaf area calculation later. 
The ambient layer in Figure 5(d) represents other objects in the image. 


(c) 


Figure 5. Background removal in top-view image: (a) original image, (b) leaf area, (c) calibration layer, and 
(d) ambient layer 
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For evaluation, we conduct some experiments with same conditions to measure the performance of 
our proposed method and compare with the results of original mahalanobis distance based method [36] for 
leaf segmentation. The comparsion results are given in Table 2. It is clear that with the growth of cucumber 
plants, the number of leaves increases. Cucumber plants with greater leaf area have a larger number of 
overlapping leaves. Consequently, the number of misclassified pixels also increases. The results in Table 2 
show that all three aspects (FNR, FPR, and ER) increase with the growth of cucumber plants. With lower 
noise points in the top-view images, all three performance indices have lower increase rates than that in the 
side-view images in all growth stages of plants (DAT). Figure 6, shows the segmentation results in side-view 
image. After removing the background of the original image in Figure 6(a), the leaf layer, calibration layer, 
and ambient layer are segmented and depicted in Figures 6(b)-(d) respectively. The leaf layer in Figure 6(b) 
shows that the number of FP and FN pixels are significantly higher than that in Figure 5(b). The results of 
our proposed MMD algorithm indicates that the averages of FNR, FPR, ER are 0.0284, 0.0201 and 0.0236, 
respectively. The average of misclassified pixels in every image is 2.36%. It reveals that our proposed MMD 
algorithm performs well and can be effectively used in most of image-processing applications. 


Table 2. The comparsion results of leaf segmentation methods 


, Original Mahalanobis distance-based method Our proposed MMD method 
DAT Image view 


FNR FPR ER FNR FPR ER 

Top-view 0.2221 0.0364 0.1071 0.0114 0.0245 0.0195 

10 Side-view 0.1402 0.0314 0.0748 0.0119 0.0155 0.0141 
Top-view 0.1890 0.0191 0.0994 0.0137 0.0144 0.0141 

ra Side-view 0.1291 0.0240 0.0698 0.0145 0.0161 0.0155 
Top-view 0.1284 0.0296 0.0730 0.0153 0.0150 0.0151 

1# Side-view 0.1775 0.0198 0.0858 0.0161 0.0161 0.0161 
Top-view 0.1628 0.0216 0.0774 0.0169 0.0167 0.0168 

16 Side-view 0.1481 0.0210 0.0773 0.0177 0.0173 0.0175 
Top-view 0.2577 0.0251 0.1303 0.0185 0.0179 0.0181 

L$ Side-view 0.1991 0.0273 0.0958 0.0194 0.0185 0.0188 
Top-view 0.1550 0.0308 0.0815 0.0202 0.0190 0.0195 

20 Side-view 0.1406 0.0240 0.0742 0.0275 0.0197 0.0228 
Top-view 0.2986 0.0402 0.1600 0.0259 0.0203 0.0225 

a Side-view 0.1779 0.0354 0.0948 0.0266 0.0208 0.0232 
Top-view 0.1049 0.0323 0.0617 0.0274 0.0214 0.0239 

24 Side-view 0.2446 0.0275 0.1264 0.0290 0.0220 0.0249 
Top-view 0.3209 0.0437 0.1793 0.0378 0.0228 0.0289 

= Side-view 0.2348 0.0376 0.1293 0.0393 0.0234 0.0299 
28 Top-view 0.2788 0.0310 0.1552 0.0455 0.0241 0.0329 
Side-view 0.1786 0.0288 0.0986 0.0517 0.0248 0.0360 
Top-view 0.2751 0.0413 0.1521 0.0637 0.0256 0.0417 

a Side-view 0.2145 0.0319 0.1147 0.0754 0.0264 0.0474 
Average 0.1990 0.0300 0.1054 0.0284 0.0201 0.0236 


(c) | (d) 
Figure 6. Background removal in side-view image: (a) original image, (b) leaf layer, (c) calibration layer, and 
ambient layer 
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According the experimental results given in Table 2, the proposed MMD algorithm achieves higher 
accuracy rate, which is 8.18% greater than that of original Mahalanobis distance based algorithm. One main 
reason of this difference is how to chose a segmentation threshold for each window of segmentation image. 
In general, the original Mahalanobis distance based methods choose a segmentation threshold for all 
segentation images, which tends to increase the segmentation error rate in homogeneous colour regions. 
While in our proposed MMD algorithm, each image being segmented is splitted into smaller regions and then 
the segmentation threshold is determined based on image local properties. As consequence, our proposed 
MMD algorithm worked well even in homogeneous colour regions. 


3.2. Leaf area estimation results 

As described above, there are two approaches for LA estimation of cucumber plants. With the top- 
view images of cucumber plants, the LA in each image is computed as given in (9). However, with side-view 
images we first estimated the number of leaves in each image and then estimated the LA based on manual 
measurements given in Table 1. In order to estimate the number of leaves in side-view images, we first find 
the skeleton of each cucumber plant and then count the nodes along these plants stems. Each node has a leaf. 
The procedure of skeleton extraction is described in Figure 7. The original image in Figure 7(a) is used for 
segmentation. Figure 7(b) shows the segmentation result of our proposed MMD algorithm. The background 
of this segmentation image is changed to white colour and given in Figure 7(c). The binary image is depicted in 
Figure 7(d). By using Zhang-Suen algorithm [48], the main stems of cucumber plants are given in Figure 7(e). 

As described in Figure 8, there is no significant difference between the estimated and measured 
length of the main cucumber stem. The relative error of these two parameters varies from 4.77% to 8.47 %, 
which indicates that our proposed algorithm can achieve high accuracy in measuring the length of the main 
cucumber stem. The highest value of relative error is only 8.47% when DAT is 16 days. 

Table 3 provides a summary of LA estimation results. It indicates that the relative error of our 
proposed approach increases from 3.76% to 13.00%. Having larger overlapped leaf regions in top-view 
images, the relative errors of leaf area estimation with the same cucumber plant are normally higher than that 
in side-view images. 


Figure 7. Schematic illustrating the process of skeleton extraction: (a) original image, (b) segmentation 
image, (c) white background segmentation image, (d) binary image, and (e) skeletion of cucumber plant 


E Estimated length of main stem(L,) 


Measured length of main stem(L,,,) 


Figure 8. Estimated lengths of main stem 
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Table 3. The comparision between top-view and side-view images prosessing for LA estimation 
DAT Image view LA. LAm RE (%) 


10 Top-view 198.6 4.20 
Side-view 199.5 20.3 3.76 
12 Top-view 242.5 5.01 
Side-view 2443 °553 431 
14 Top-view 352.6 371.1 4.99 
Side-view 356.2 i 4.02 
16 Top-view 443.0 4.30 
Side-view 440.6 4623 4.82 
18 Top-view 502.6 6.02 
Side-view 498.6 DLR 6.77 
20 Top-view 577.6 613.5 5.85 
Side-view 565.3 : 7.86 
22 Top-view 712.5 6.61 
Side-view 7096 72? — 699 
24 Top-view 850.3 6.30 
Side-view 845.3 ie 6.85 
26 Top-view 1028.3 7.71 
Side-view 10459 11142 613 
28 Top-view 1210.7 8.85 
Side-view 12075 1783 9.09 
30 Top-view 1372.9 13.00 
Side-view 1382.7 13780 12.38 
Average 6.63 


4. CONCLUSION 

This paper presents two new efficient approaches for cucumber leaf area estimation. In the first 
approach, the LA of cucumber plants in top-view images was simply estimated based on the calibration area. 
In this way, our proposed approach achieves a high accuracy rate while mitigating the manual camera 
calibration. In the second approach, three processes including background eliminating, skeleton extraction, 
and leaf number estimation are sequentially applied for cucumber LA estimation. Our experimental results 
demonstrate the effectiveness of our proposed approaches. Limitation of the paper is the fact that the quality 
of captured images be influenced by ambient light noise. The appearance of some overlapped leaf areas is 
also the main reason of high relative error. This work, however does not consider the real time estimation of 
leaf area. Therefore, our future plan is to enhance these issues with some realistic LA estimation algorithms. 
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